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Abstract 



In this letter we focus on the calculation of f^, , the amplitude of the (cc) — > 
(gluons) — ► r\ transition, which provides the magnitude of the contribution 
of the Cabbibo favored b — * ccs elementary process to B — > rj'K decay. It 

(c) 

is found that f\ J = —12.3 ~ —18AMeV on the scale of m c . This number 
is in strong contradiction with the estimations of Halperin, Zhitnitsky and 
Shuryak, Zhitnitsky but almost in agreement with phenomenological analysis 
of Feldman et al., Petrov and the estimations of Ali et al. 
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I. INTRODUCTION 



Recently there is a great theoretical interest P|-|i"3"| on the recently released experimental 
data on the branching ratios of the decays of B — > Krf I0J, |J: 

B(B ± -> rj'K^) = (6.5+};| ± 0.9) x 10" 5 , (1.1) 
B(B° -> 7/ AT ) = (4.7j£o ± 0.9) x 10" 5 . (1.2) 

In the Standard Model the Cabbibo favored b — > ccs elementary process may be followed by 
conversion of cc pair into rf through the gluons. The amplitude of this process is described 
by 

M = %W>1 < V'(p)\ci^c\0 X K(q)\s lfl b\B(p + q) > . (1.3) 

Here Gf is the weak coupling constant, V c b, V* s are Kobayashi-Maskawa matrix elements, 
a\ = 0.25 phenomenological number obtained by a fit (see [§] for the references). The matrix 
element 

<v'{p)\n^ 5 c\0>=-if^ Pfl (1.4) 

is non-zero due to the virtual cc — > gluons transitions. Certainly, this matrix element is 
suppressed by the l/m 2 c factor. However, due to strong nonperturbative gluon fields together 
with the Cabbibo favored b —>■ c transition the suggested mechanism (|1.3|) can be expected 
to compete appreciably with the other mechanisms of the B — > Kr]' process P|L2fl. 

If we assume the dominance of the mechanism (|1.3|) the branching ratio is written in 
terms of fjy as 

/ Ac) \ 2 

Br[B Kr() ~ 3.92 ■ 10" 3 ■ . (1.5) 

Using the data Q it is found fjf ~ 140 MeV ("exp"). 

This value perfectly coincides with estimation of Halperin and Zhitnitsky [[| : 

= (50 - 180) MeV. (1.6) 

On the other hand, a recent phenomenological study placed a bound on , namely 

-65 MeV < f { f < 15 MeV, with ffi being consistent with zero by analyzing the Q 
evolution of the r] ,r y form factor |7j], and more recently it was estimated from observed ratio 
of J/ip decay to v( and n c the value of fffi = — (6.3 ± 0.Q)MeV p. Other similar estimation 

which leads to < 12MeV was made in ||13|| . Ali et al. considered the complete am- 

plitude for the exclusive B— meson decays, including n'K channels, where it was combined 
the contribution from the process b — > s(cc) — * s(gluons) — > sr/( > with all the others arising 
from the four-quark and chromomagnetic operators, as detailed in their papers ||, @. Their 
estimations gave \ f^\ — 5.8 MeV f| and fffi = —3.1 (—2.3) MeV (for m c in the range 1.3 
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-1.5 GeV) || in agreement with the analysis j7|. They stressed the importance of the sign 

(c) 

of /V and found a theoretical branching ratio in the range 

B(B — > r/K) = (2-4) x 1(T 5 , 

which is somewhat smaller than the experimental one Ql.ip . The similar analysis made in 
|ICfl leads the authors to conclude that fffl = —50MeV may provide the explanation of the 
data. 

Having this situation, it is important to recalculate f!$ to clarify the mechanism of 
B — > Krf decay in the similar framework performed by Halperin and Zhitnitsky ||. 



II. CALCULATIONS OF 

The symmetry of the classical lagrangian may be destroyed by quantum fluctuations 
nj-(L6|. In gauge theories the axial anomaly arises from noninvariance of the fermionic 



measure against axial transformations in the path integrals of the theory |T7j (see also ref. 
concerning higher-loop corrections). The present problem is intimately related with 



this phenomena. 

In the following we will work only with the Euclidian QCD[|. In the Euclidean QCD the 
axial anomaly in the light quark axial current in chiral limit reads 

d^fm^f = -*yf^G, (2.1) 

where ipf is the light quark field (/ = u,d,s) and g the QCD coupling constant. 2GG = 
e nv\crQa^Qa^^ w here G a is the gluon field strength operator with a being the color index. 

The situation with heavy quarks is very different, since we must take into account the 
contribution of the mass term. The divergence of the axial current of charmed quarks has a 
form: 

a 2 

d^y^c = -ij^GG + 2mJ l5 c, (2.2) 

The first term in ( j2.2p again comes from noninvariance of the fermionic measure (or in 
other words - from Pauli-Villars regularization). The main problem here is to calculate the 



contribution from the second term in (|2.2j ). It is clear that this one is reduced to the problem 
of the calculation of the vacuum expectation value of the operator 2m c c^y^c in the presence 
of a gluon fields. 

In the path integral approach the calculation of the contribution of this term to any 
matrix element over light hadrons may be considered in sequence of the integrations. Firstly 



we are using here the convention of the Euclidian QCD: 

ix M = XEA, X M i = XEi, A M = iA E 4, A M i = ~A E i, V'A/ = TpE, i^M = V4> 7M0 = 1E4 

iMi = ilEi, 7A/5 = lEb- We will omit index E. 
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the integration over c-quark is performed, and the next step is the calculation of the integral 
over gauge gluon field and finally integral over light quarks. 

We consider here the first step - the integration over c-quarks. We define: 

< 2m c c\x)^i^c{x) >= / DcDcHm c c\x)^c(x) exp( / c f (iV + im c )c). (2.3) 



Here iV = "i^iid^ + gA p ) and we introduce the operator P p and p M which are defined in 
the coordinate space as < x\P p \y > = iV p 5(x — y) and < x\p^\y > = id p 5(x — y). It is clear 
that the formal answer for the path integral ( (2.3|) can be written in the form: 



< 2m c c^(x)75c(x) >= 2m c det\ \P + im c \ | < x\Tr^§{P + im c ) l \x > (2.4) 
det\\P + im c \ \ must be regularized in the standard manner as 

, ii #, • i' i n (P + im c )(p + iM) n 

det \\P + im c \\ -> det , 

(p + im c )(P + iM) 

where M is the regulator mass. Eq. (|2.4j ) must be a gauge invariant function of the gauge 
field A and therefore must be expressed through the gluon field strength tensor and their 
covariant derivatives. 

We will follow the operator method proposed by Vainshtein et al in the same line 



as in [[|. The key ingredient of this method is based on an assumption of a possibility of 
an expansion of ( |2.4| ) over We will take into account 0{g 2 G 2 ) and 0(g 3 G 3 ) terms in the 
calculations of ( |2.4| ). We start from the calculation of 

H(x) = 2m c < x\Tr^ 5 (P + im c )~~ l \x > = (2.5) 

2 



-2im\ < x\T n5 (P 2 + m 2 c + ^agG)' 1 ^ > = H 2 {x) + H 3 {x) + 0{g*G 4 ), (2.6) 



where 

H 2 (x) = -2%m\ < x\T n5 (P 2 + m 2 c )-%oG{P 2 + ///f) v9 -aC > (2.7) 

and 

H 3 (x) = 2im 2 <x\Tr l5 (P 2 + ml)-' z ^aG(P 2 + m 2 c )- L ^aG(P' z + m 2 c )-^aG]\x> . (2.8) 



Here aG = a^G^, a^ v = f [7^,7^- 

It is straightforward to calculate H 3 (x), since we may neglect the noncommutativity of 
the operators in (|2.8| ) and replace P operator by p in (|2.8|). In that case we may use the 
evident formulas, 

< x \ 1 \ x >= f d \ 1 = ( 2 V(n - l)(n - 2)m 2 (™- 2 ))- 1 , 
> 2 + m 2 )™' 7 (2tt) 4 (p 2 + m 2 ) n ' ' 

Tr 75 (aG) 3 = i2 5 tr c GGG, 
where GGG = G fJiU G up G pp . We get 
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H *i x ) = - kJ] 2 foj bc G a &G c (2.9) 

c 



2 4 3vr-m 



The calculation of H 2 (x) needs much more efforts. First of all we represent the aG(P 2 + 
m 2 ) -1 in the form 



CJ 



aG(P 2 + m 2 c )' 1 = {P 2 + m 2 c )- l aG + (P 2 + m 2 c )~ 2 [P\ aG] + (P 2 + m 2 )" 3 [P 2 , [P 2 , cxG]] + ... 

(2.10) 

The routine calculations of the commutators in ( p.lOj ) lead to 

[P 2 ,aG] = V 2 aG + 2iP a V a aG, (2.11) 

[P 2 ,[P 2 ,aG]] = VVG + 2zP Q (V Q V 2 + V 2 V a )aG + {2i) 2 P G Pa V a oG + {2i) 2 P a P p V P V a <rG. 

(2.12) 



Other higher commutators lead to the terms order 0(G 3 ) in the expansion Q2.10D and may 



be neglected. Following the arguments of [19] we may neglect also the terms which contains 



single operator P M . The reason is that the matrix elements 

<x\(P 2 + m 2 c y n P ll \x>~ V^G 2 . 

By using the Bianchi identity it is easy to show that 

V 2 G^ = i(G au G afM + G a ^G ua ) — V v V a G aAt — V M V a G ua (2-13) 

It is evident that the term V 4 crG is order of 0(G 3 ) and may be neglected. Collecting all of 
the 0(g 2 G 2 ) and 0(g 3 G 3 ) terms in H 2 (x) in (|277]), we get 

H 2 (x) = ^G>& + i^^U^G^ (2.14) 



Finally 



■ 2 • 3 

H(x) = ^—G a G a - J abcG a G b G c (2.15) 

v > 2 4 tt 2 2 5 37T 2 my v ; 



We neglect here the small contributions of the terms like V^V a G ua . 

As expected, the first term in H(x) cancels with the first term in ( |2.2|) , which is the 
contribution from noninvariance of the measure and the rest part leads to the divergence of 
the c-quark axial current in the form 

• 3 

< a M ct(x) 7/i7feC (x) > = - n ,l 9 M fa bc G a G b G c . (2.16) 

c 



2 5 37r i, m 



We would like to stress an attention that our answer for < d^c^x)^^ jsc(x) > is 6 times less 
than was calculated by Halperin and Zhitnitsky ||. 
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We apply this result to the calculation of the f^f . The analogous quantity f^., , which 
is defined in the similar way as f$ in( |1.4| ), is 

mlf^ =<0\^- 2 G a G a \ V '>. (2.17) 



By definition (|1.4|) must be calculated from: 



m\f^ = < Ol ^^ fabcG^G^' > (2.18) 

The phenomenological way of the estimation of the fty is the application of the QCD+QED 
axial anomaly equation together with data on rf — > 27 decay leads to 

/$° = Q3.QMeV, (2.19) 



which was used in ||. In the DP chiral quark model |2T|, which was successfully checked 
by the calculation of the axial anomaly low-energy theorem in chiral limit, the matrix 
element in ( |2.18| ) [ 22" ] , [ 2D | may be reduced to the calculation of the matrix element in 17| ) 
with an additional factor |22|, |2pJ. Here p is the average size of the QCD vacuum 

instantons. Phenomenological analysis, variational and lattice calculations showed that 

p = 1/3 fm (2.20) 

So, the ratio of Eqs. (|2.18 ) and ( 2.17 ) is equal in this model to: 

4t = = = 0.1 2.21 

5p 2 Qm 2 c v ' 

By taking into account the estimation ( [2.19|) (we use m c (pi ~ m c ) ~ 1.25 GeV on the scale 
Pi ~ m c for the numerical estimates), we find 

ff = — QMeV (2.22) 

This number is close to the one of [§, \fft\ = 5.8MeV and the sign and the order of the 
value coincide with the estimations of ||§,|13|]. 

Recently, Shuryak and Zhitnitsky M performed direct numerical evaluations of the var- 
ious correlators of the operators g 2 G a G a , g 3 f a bcG a G b G c in the Interacting Instanton Liquid 
Model(IILM). Their calculations lead to: 

< 0\g 2 G a G a \ V ' > = 7GeV 3 (2.23) 
(which leads to ffl = 48.3 MeV) and 

\ m t \ 1 

(1.5 ~ 2.2)GeV 2 (2.24) 



1 < 0\g 3 f abc G a G b G c 


rf 


> 




< 0\g 2 G a G a \r]' 


> 
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The later is somewhat large than their simple estimate for this ratio of matrix elements 

12 1 



< — >^ (1 ~ 1.5)GeV 2 (2.25) 



5 p 



They concluded that 



1$ 



(c)| 



(«)l 



1.47 ~ 2.11 (2.26) 



On the other hand, with the use of ( |2.24j ) and the abovementioned factor —1/6 in ( |2.16 ) we 
arrive at 

Ac) 

12L = -0.17 0.25. (2.27) 

This ratio gives = —8.2 ~ —12.3 MeV at the scale of the size of the instanton /i 2 ~ /o -1 . 
The abovementioned experimental numbers ( |1 . 1| ) are given at the scale fii ~ m c , which 
is different from the scale of this instanton calculation. The account of the anomalous 
dimension of the g 3 GGG operator leads to correction 

4, c Vi ^ rn c ) ^ 1.5/^(^2 ^ p- 1 ), (2.28) 
Q. The account of this scale factor leads to 

- m c) = -12.3 ~ -18.4Me^ (2.29) 



Hence, using (|2.24j) , the result of more sophisticated calculations of Shuryak and Zhitnitsky, 



we get the number ( |2.29| ) which is 2-3 times larger than simple estimation ([2.22 ). 



These numbers fl2.22|) , fl2.29|) are in agreement with the phenomenological bounds |7|.|T3 



and almost in agreement in the sign and the value with but six-ten times less than the 
estimations given by @(see (|1.6|) ) and also [§[] (see ( |2. 26| ) ) . 



By using the numerical analysis of the branching ratio for —>■ rj'K given at 
(Fig. 17 of ) we expect that the value of given in fl2.29p may provide a more satis- 



factory explanation of the experimental data|. We reserve this investigation for the future 
publication. 
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quite recently Cheng and Tseng [11| concluded impressive good explanation of these data using 
our value (B) for &\ 
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